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Abstract. We describe observations of tropical stratospheric water vapor q that 
show clear evidence of large-scale upward advection of the signal from annual 
fluctuations in the effective “entry mixing ratio” q E of air entering the tropical 
stratosphere. In other words, air is “marked,” on emergence above the highest 
cloud tops, like a signal recorded on an upward moving magnetic tape. We 
define q E as the mean water vapor mixing ratio, at the tropical tropopause, of 
air that will subsequently rise and enter the stratospheric “overworld” at about 
400 K. The observations show a systematic phase lag, increasing with altitude, 
between the annual cycle in q E and the annual cycle in q at higher altitudes. The 
observed phase lag agrees with the phase lag calculated assuming advection by the 
transformed Eulerian-mean vertical velocity of a q E crudely estimated from 100-hPa 
temperatures, which we use as a convenient proxy for tropopause temperatures. 
The phase agreement confirms the overall robustness of the calculation and strongly 
supports the tape recorder hypothesis. Establishing a quantitative link between 
( 1e and observed tropopause temperatures, however, proves difficult because the 
process of marking the tape depends subtly on both small- and large-scale processes. 
The tape speed, or large-scale upward advection speed, has a substantial annual 
variation and a smaller variation due to the quasi-biennial oscillation, which delays 
or accelerates the arrival of the signal by a month or two in the middle stratosphere. 
As the tape moves upward, the signal is attenuated with an e-folding time of about 
7 to 9 months between 100 and 50 hPa and about 15 to 18 months between 50 and 
20 hPa, constraining possible orders of magnitude both of vertical diffusion K z and 
of rates of mixing in from the extratropics. For instance, if there were no mixing 
in, then Ii z would be in the range 0.03-0.09 m 2 s _1 ; this is an upper bound on K z . 
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1. Introduction 

It is well known that the annual cycle in water va- 
por mixing ratio q at the tropical tropopause appears 
to spread sideways to middle latitudes, with a certain 
phase lag. This is shown, for instance, by data from 
radiosondes [Mastenbrook and Oilmans, 1983; Hyson , 
1983], from the Stratospheric Aerosol and Gas Exper- 
iment (SAGE) II [McCormick et al, 1993], and from 
aircraft observations [Hintsa et ai, 1994]. For the trop- 
ical lower stratosphere itself, Mastenbrook and Oilmans 
[1983] called attention to Mastenbrook? s [1968] sound- 
ings at Trinidad (11°N, 60°W), which showed a deepen- 
ing dry layer during northern winter and spring. Kley 
et ai [1979] noted a decrease of q with height from 
the tropopause to 19 km in a balloon sounding above 
Brazil in August. They suggested a few possible expla- 
nations, one of which attributed the minimum at 19 km 
to upward advection of the low q values that presum- 
ably existed at the tropopause 6 months earlier. Other 
in situ observations of q have shown layered structures 
in the tropical lower stratosphere at various times [ Kley 
el ai, 1982; Kelly et ai, 1993; Vomel et ai, 1995; Wein- 
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stock et ai f 1995; J. Ovarlez et al., In situ water vapour 
measurement, and a case study of the drying mecha- 
nism of the tropical lower stratosphere, submitted to 
Quarterly Journal of the Royal Meteorological Society , 
1995, hereinafter referred to as Ovarlez et al., submitted 
manuscript, 1995]. 

At a NATO workshop on stratosphere-troposphere 
exchange in September 1993, one of us (M.E.M.) ar- 
gued that such tropical lower stratospheric q features 
should be part of a larger spatio-temporal pattern with 
a very simple phase lag structure. Air passing the tropo- 
pause, and carried upward in the large-scale strato- 
spheric mean circulation, should be marked with the 
entry mixing ratio </e in the same way as a magnetic 
tape is marked by the recording head. The variable- 
speed “motor” driving the “tape'’ is the well known 
extratropical wave driving of the large-scale mean cir- 
culation, in which the extratropical stratosphere and 
mesosphere act together as a global-scale suction pump 
on the tropical lower stratosphere [e.g., Holton et ai, 
1995, and references therein], with seasonally and inter- 
annually variable strength [Rosenlof and Holton, 1993; 
Yulaeva et ai, 1994]. 

Direct evidence in support of this tape recorder hy- 
pothesis has recently become available from observa- 
tions of q by the Microwave Limb Sounder (MLS) in- 
strument on the Upper Atmosphere Research Satellite 
(UARS) [Mote ct ai, 1995]. The larger-scale view re- 
vealed by the MLS data showed what looked like the 
signature of upward advection of q features by the 
mean circulation, seemingly little degraded by mix- 
ing. Mote et al. [1995] checked this by using observed 
ozone, water vapor, and temperature, together with a 
state-of-the-art radiation scheme, to calculate a trans- 
formed Eulerian-mcan (TEM) circulation and running 
two-dimensional (2-D) backward trajectories from var- 
ious latitudes and heights in the stratosphere down to 
100 hPa. The q for the entire trajectory was then es- 
timated as a saturation value q s based on gridded 100- 
hPa temperatures. The q values observed by MLS var- 
ied with time in a way sufficiently like that of the 2-D 
trajectory calculation to confirm, or at least to demon- 
strate consistency with, a tape-recorder-like connection 
between those q values and the annual cycle in tropo- 
pause temperature. 

There are a number of questions about the tape 
recorder hypothesis that were not answered by Mote et 
ai [1995]. They found a large discrepancy in the am- 
plitude of the observed q variations and the estimated 
tape signal. What was the cause of this amplitude dis- 
crepancy? Would observations at altitudes below 46 
hPa, beyond the range of the standard MLS retrieval 
algorithm, show the same kind of pattern? Can mixing 
effects be quantified? What is the effect of the quasi- 
biennial oscillation (QBO)? 

In this paper we address these questions with the 
help of additional data and calculations. Some of the 
data come from three instruments aboard UARS: MLS, 
with an improved water vapor retrieval, the Halogen 
Occultation Experiment (IIALOE), and the Cryogenic 


Limb Array Etalon Spectrometer (CLAES). Also used 
are data from the Stratospheric Aerosol and Gas Exper- 
iment (SAGE II) aboard the Earth Radiation Budget 
Satellite, from the NOAA Lyman a hygrometer aboard 
the NASA ER-2 aircraft, and from radiosondes. The 
sources of data are described in section 2. In section 
3 we use these data and some calculations to exam- 
ine in more detail whether the tape recorder hypothesis 
does indeed explain the observations. In particular, we 
look at the QBO effects (section 3.1), at some possi- 
ble causes of the amplitude discrepancy (section 3.2), 
and at the role of mixing in attenuating the tape signal 
(section 3.3). Section 3.2 focuses on the relationship 
between 100-hPa temperatures (as obtained, for conve- 
nience, from the United Kingdom Meteorological Office 
(UKMO)) and observed q in the stratosphere. For this 
purpose we define q e ns the mean water vapor mixing 
ratio, at the tropical tropopause, of air that will subse- 
quently rise and enter the stratospheric ‘ overworld at 
around 400 K or 83 hPa. Air that reaches the tropo- 
pause but fails to get into the overworld, by returning 
to the troposphere either downward or poleward, for 
example, will therefore not affect q^. 

2. Data Sources and Treatment 

Because each instrument has different observational 
characteristics, we have used different methods of pre- 
senting the data from each instrument. The choice of 
latitude limits to define the tropics is not very criti- 
cal. As was shown by Carr et ai [1995], q has small 
meridional variations between about 15°S and 15° N, so 
temporal variations tend to be in phase. Poleward of 
those approximate limits, q has a sharp gradient beyond 
which annual variations are small and are not as well 
correlated with tropical annual variations. 

Gaps in the UARS data record appear at various 
points. The longest gap occurred during the period 
June 2 to July 15, 1992, when problems with the drive 
for UARS’s solar array kept the instruments from func- 
tioning. The MLS water vapor radiometer was turned 
off for most of this period, and IIALOE was turned off 
for the first few weeks. 

Because of methane oxidation, which produces wa- 
ter, there is generally an increase of water vapor with 
height in the stratosphere, and this increase makes it 
difficult to relate annual variations at different lev- 
els. This difficulty can be alleviated by using twice 
methane plus water, 2(CH4)+(H 2 0), in place of water, 
because methane oxidation has a conversion efficiency 
very close to 2 [e.g., Dessler et ai 1994]. The quan- 
tity 2(CH 4 )-|-(H 2 0) (hereinafter called H , as suggested 
by Jones et ai [1986]) is the variable parr of “total 
hydrogen,” 2(CH 4 )+(H 2 0)+(H 2 ) 5 since (H 2 ) is nearly 
constant at about \ part per million by volume (ppmv) 
in the lower and middle stratosphere [e.g., Dessler et 
ai 1994]. HALOE and CLAES both retrieve methane 
and we can plot H for those instruments (though be- 
cause of problems with CLAES H 2 0 we do not do so 
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for CLAES). MLS and SAGE II do not, however, re- 
trieve methane, so we must follow a somewhat cruder 
approach and subtract the time mean vertical profile, 
where the time mean is taken over the entire record in 
each case. 

2.1. Microwave Limb Sounder (MLS) 

MLS vertically scans the Earth’s limb and measures 
atmospheric thermal radiance at millimeter wavelengths, 
from which vertical profiles of atmospheric parameters 
are retrieved. The measurement technique is described 
by Waters [1993], and the instrument is described by 
Barath et ai [1993]. Stratospheric water vapor pro- 
files were obtained from the MLS 183-GHz radiome- 
ter, which operated from late September 1991 through 
late April 1993. The tropics were observed continuously 
(with few exceptions) during this period, and approx- 
imately 240 profiles between 12°S and 12°N were ob- 
tained each day. The estimated accuracy is better than 
20%, and individual profile points have precisions of 
0.2 to 0.4 ppmv [Lahoz et al, 1996]. Since hundreds of 
profiles are averaged for our analyses, the observational 
results presented here should have very good precision. 

For the standard version 3 retrieval, the vertical pro- 
files are piecewise linear, with nodes (grid points) sepa- 
rated by a factor of 10 1 / 3 in atmospheric pressure. This 
gives levels at 100, 46, 22, and 10 hPa in the lower 
stratosphere, separated by approximately 5.3 km. A cli- 
matology based on Limb Infrared Monitor of the Strato- 
sphere (LIMS) data [Gille and Russell , 1984] is used as 
a priori data for the retrieval, and the a priori stan- 
dard deviation is 2 ppmv. A strict opacity criterion is 
imposed that prevents radiances being used if the atmo- 
sphere is too opaque at 183 GHz along the line of sight; 
this opacity criterion must be imposed because the ver- 
sion 3 retrieval is linear, that is, radiances are presumed 
to be proportional to the concentration of water vapor 
in the tangent layer. As a result of the opacity crite- 
rion, the 46-hPa level is the lowest level for which water 
vapor is retrieved. 

To overcome this limitation, we reprocessed tropical 
MLS data for a number of days (at approximately 10- 
day intervals) using a different, nonlinear procedure. 
We used a grid with twice as many levels, that is, with 
nodes separated by a factor of 10 1 / 6 in atmospheric 
pressure, or about 2.7 km; this grid is more closely 
spaced than the theoretical resolution of MLS, which 
is about 3.5 km, but we can still retrieve a profile (as 
described below). The finer grid has levels at 100, 68, 
46, 32, 22, 15, and 10 hPa in the lower stratosphere. As 
a priori data for the retrieval we used a revised climatol- 
ogy based on SAGE II, and the a priori standard devia- 
tion was 3.5 ppmv. Lest the monthly mean climatology 
introduce an artificial seasonality, we used the January 
climatology as the a priori throughout the year. Instead 
of imposing the opacity criterion, we used a nonlinear, 
iterative retrieval that (fortunately) requires only a few 
iterations in the tropical lower stratosphere. At each 
iteration we recalculated only a few of the processing 
coefficients instead of the full set. We ascertained that 


this is a reasonable approximation to the full nonlinear 
retrieval by comparing the two approaches for one or- 
bit; the agreement was excellent in the tropical lower 
stratosphere, although it was poor in the polar regions 
and in the mesosphere. The a priori covariance matrix 
has some nonzero off-diagonal elements that smooth the 
retrieved profile just enough to prevent computational 
instabilities. (An element Sij of the covariance matrix 
represents the covariance between errors at retrieval lev- 
els i and j ; use of a diagonal covariance matrix would 
mean that errors at different levels are assumed to be 
uncorrelated. With a grid finer than the instrument’s 
theoretical vertical resolution, no unique solution exists, 
but if off-diagonal elements of the covariance matrix are 
nonzero, then the retrieval is lightly constrained to pro- 
duce a smooth profile rather than a jagged profile.) 

The results of this procedure should be regarded as 
preliminary for at least two reasons. First, they use 
MLS version 3 temperature and pressure values, which 
have known minor problems (these will be reduced in 
future versions of the data). Second, the detailed error 
analyses that have been carried out on version 3 data 
have not been performed on the data shown in this pa- 
per. 

We are confident, however, that the results given here 
represent a real improvement over version 3. The error 
ratio, that is, the ratio of retrieved error to a priori stan- 
dard deviation, is small even at 68 hPa and usually also 
at 100 hPa, indicating good data (and minimal reliance 
on the a priori data) at those levels. The finer verti- 
cal resolution is crucial to resolving the tape signal, as 
was apparent from the weak visual connection between 
extrema at 46 and 22 hPa in Figure la of Mote et ai 
[1995]. 

2.2. Halogen Occultation Experiment (HALOE) 

HALOE (described by Russell et al. [1993]) is a solar 
occultation instrument that, because of the inclination 
of the UARS orbit, samples the atmosphere about 30 
times per day, half at one latitude at sunrise and half at 
another latitude at sunset. These two observation lat- 
itudes change with time as the UARS orbit precesses, 
so the tropical observations shown are actually taken 
at irregular intervals (see Figure 2) but are binned by 
month for Plate 1. H ALOE’s orbital characteristics pre- 
vent it from viewing the tropics at all in some months. 
But HALOE has an important advantage in that it does 
not require an absolute radiance calibration in order to 
make its measurements. Instead, it calculates the ratio 
of the solar intensity measured through the atmosphere 
at the Earth’s limb to the unattenuated solar intensity 
measured above the atmosphere. Consequently, the ex- 
periment provides very precise data that are well suited 
for time trend studies. We use HALOE version 17 data; 
the version 17 retrieval algorithm represents a marked 
improvement over earlier versions. The effective verti- 
cal weighting of HALOE measurements is about 2 km, 
but the data are shown on a vertical grid of 1.3 km for 
Plate lb and of 0.2 km for Figure 1. 
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2.3. Cryogenic Limb Array Etalon Spectrometer 
(CLAES) 

CLAES (described by Roche et ai [1993]) is a limb- 
viewing spectrometer that measures the infrared emis- 
sion by several gases and by aerosols. In order to reduce 
overlap and interference between gases, CLAES uses 
high spectral resolution, which requires that the instru- 
ment be cooled to temperatures of 15-130K (by solid 
neon and solid carbon dioxide) to allow narrow-band 
measurements to be made with high signal-to-noise ra- 
tios. CLAES provides data every 500 km along the 
orbit, with a vertical resolution of 2.5 km, spanning 50 
km in altitude from a nominal lower level at 10-15 km. 
Because of the opacity of atmospheric features at low 
altitudes and the effects of the Pinatubo aerosols, the 
retrievals at altitudes below 46 hPa are not as reliable 
as those at higher altitudes. We use CLAES version 
7 water vapor data, omitting values at 46 hPa before 
July 1992 because they appear to be contaminated by 
aerosol. 

2.4. Stratospheric Aerosol and Gas Experiment 
(SAGE II) 

The characteristics of the SAGE II water vapor re- 
trieval and measurements are described by Rind et ai 
[1993]. Like HALOE, SAGE II is a self-calibrating solar 
occultation instrument. It has a nominal vertical res- 
olution of 1 km and provides the most extensive avail- 
able water vapor data set, extending in time from 1985 
to 1991, in latitude from about 70°N to 70°S and in 
altitude from the surface to 120 km. SAGE II uses 
the 0.94 /im wavelength near-infrared channel to de- 
rive water vapor mixing ratios. Because many aerosol 
particles are a substantial fraction of a wavelength in di- 
ameter and therefore strongly scatter the near-infrared 
radiation, too much aerosol makes the retrieval of wa- 
ter vapor impossible. In the tropical lower stratosphere 
there is always enough aerosol between 20 and 25 km to 
make retrievals very noisy at best. Following the erup- 
tion of Mount Pinatubo in June 1991, the enhanced vol- 
canic aerosol catastrophically interfered with the SAGE 
II water vapor retrievals, which were therefore discon- 
tinued. Note incidentally that the aerosol problem has 
little effect on the midlatitude results of McCormick et 
ai [1993], which were for the years 1986-1988. 

The SAGE II data set available to us extended from 
January 1986 to May 1991, 65 months of data. Trop- 
ical water vapor data in the lower stratosphere show 
a steady decrease that is not fully understood (E.- 
W. Chiou, personal communication, 1994). As with 
HALOE, SAGE II does not sample the full range of 
latitudes in every month and, in particular, it samples 
the tropics for only a few days each month. 

In attempting to increase the signal-to-noise ratio 
enough to find the tape signal, we took advantage of 
the length of the SAGE II record and averaged the years 
together as follows. After discarding water vapor obser- 
vations whose retrieved error exceeded half the value of 
the water vapor mixing ratio itself, we binned the ob- 


servations by latitude in 5° bins, detrended the data for 
each altitude grid point in each bin, and then binned the 
observations by calendar month to find a mean annual 
cycle. Then we calculated the annual and semiannual 
harmonics and added them to produce a smoothed an- 
nual cycle. Then those bins and altitudes for which the 
original mean annual cycle is well represented by the 
smoothed annual cycle were used to form a tropical av- 
erage from 15°S to 15°N. “Well represented” here means 
that the correlation coefficient between the original an- 
nual cycle and the smoothed annual cycle exceeded 0.5. 
This excludes most of the bins in the aerosol layer, that 
is, between about 20 and 25 km. Note that in regions 
where a well-represented annual cycle does not exist, 
this method will produce a temporally constant value 
of zero. 

2.5. Lyman a Instrument Aboard the ER-2 

In October 1994 the NASA ER-2 aircraft made a few 
flights near the equator, as part of the Airborne South- 
ern Hemisphere Ozone Experiment/Measurements for 
Assessing the Effects of Stratospheric Aircraft (ASHOE- 
/MAESA) campaign. The NO A A Lyman a hygrometer 
[Kelly et ai, 1989] measures total water by first vaporiz- 
ing any ice, then photodissociating the vapor, and then 
finding the intensity of the resulting fluorescence. In 
the absence of cloud ice, the hygrometer simply mea- 
sures water vapor. Because no noticeable cloud ice was 
present at the altitudes examined here, the results we 
present show water vapor only. 

3. Results and Discussion 

Plate 1 shows time-height sections of deviation from 
time mean water vapor mixing ratio q for MLS (Plate 
la) and for SAGE II (Plate lc), and of H= 2(CH 4 ) + 
(H 2 O) for HALOE (Plate lb). For MLS we show daily 
averages at roughly 10-day intervals for the new, non- 
linear retrieval, from October 1991 to April 1993. For 
HALOE we show monthly means from October 1991 
to July 1995. The SAGE figure shows two cycles of 
the annual and semiannual harmonics of the 1986-1991 
detrended mean q , as was discussed in section 2.4. 

The new MLS q (Plate la) clearly suggests the up- 
ward advection of successive minima and maxima — far 
more clearly than the standard version 3 — and indicates 
that the transit time from 100 hPa to 46 hPa is 6 months 
or more. (Recall that altitudes below 46 hPa are be- 
yond the range of the version 3 retrieval.) With the 
new retrieval, MLS not only provides information on 
the lowest 5 km of the stratosphere, but also provides 
higher vertical resolution throughout the stratosphere. 

HALOE has twice the vertical resolution of even the 
new MLS retrieval and retrieves methane as well, hence 
H . As discussed above, H shows the tape signal di- 
rectly, without having to subtract the time mean. Plate 
lb shows minima and maxima in H traveling upward as 
in Plate la, and confirms the long transit time from 100 
hPa to 46 hPa. Plates la and lb also suggest that q fea- 
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Plate 1 . Time-height sections of (a) MLS water vapor mixing ratio q shown as the deviation from 
the time-mean profile, between 12°S and 12°N; (b) HALOE H =2(CH4)-f(H20), the variable part 
of total hydrogen, between 12°S and 12° N; (c) SAGE II water vapor mixing ratio q , annual and 
semiannual Fourier harmonics for January 1986 to May 1991, between 15°S and 15°N (retrieval 
affected by aerosol layer, 20-25 km); and (d) TEM 2-D back-trajectory calculation of q (see 
text). For each panel, tick marks on the ordinate indicate the vertical grid points used, and 
the color scheme is as indicated in the color bar to the right. This figure is available on the 
Internet at web site http://www.damtp.cam.ac.uk/atmos-dynamics, or by anonymous ftp from 
ftp.damtp.cam.ac.uk, cd pub/papers/mem, get tapel.ps. 
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Figure 1. Vertical profiles of water vapor q measured 
by the NOAA Lyman a hygrometer flying on the ER-2 
aircraft (diamonds) and by HALOE (solid line). The 
Lyman-a profile is from a tropical dive of the aircraft 
on October 29, 1994, in the latitude range 2.2°N to 
2.2°S and longitude range 158. 0°W to 158. 6°W. The 
HALOE profile is an average of measurements made 
between November 1, 1994, and November 13, 1994, 
in the latitude range 5°N to 5°S and longitude range 
180° to 140°W. 


tures travel upward relatively quickly during northern 
winter (about 0.4 x 10 -3 m s -1 ) and relatively slowly 
during the rest of the year (about 0.2 x 10“ 3 m s -1 ). 
We will return to this point in section 3.1. The ampli- 
tude of the tape signal diminishes sharply between the 
tropopause and about 50 hPa but diminishes much less 
between 50 and 20 hPa. The high and low H values 
seem to be rather variable from month to month, but 
this may be due to poor sampling; see Figure 2 below. 
We will draw new inferences about mixing from these 
observations in section 3.3. 

In SAGE II (Plate lc), the tape signal is not as clear 
as in MLS and HALOE, because of the difficulty of ob- 
serving q in the aerosol layer (20-25 km or so). Rela- 
tively dry air leaves the tropopause in a period centered 
around February, and relatively moist air leaves in a pe- 
riod centered around August; some months later these 
features are still evident at 19 km. The tape signal, 
however, is almost invisible in the aerosol layer. Only 
a vestige seems to be visible as a dry (moist) anomaly 
at 22 km in December (June), consistent with MLS, 
HALOE, and the transformed Eulerian mean (TEM) 
calculation. Above 25 km the dry anomalies in May and 
the moist anomalies in November compare well with the 
corresponding features in Plates la and lb, consistent 
with the expectation that SAGE II should see the tape 
signal re-emerging above the aerosol layer. 

We also show (Plate Id) results from the 2-D trajec- 
tory calculation described by Mote et al. [1995] and 
briefly summarized here. Monthly mean TEM stream 
functions are calculated for 1992 and 1993 as described 
by Rosenlof [1995]. The method uses radiative heat- 
ing rates based on observed temperatures, ozone, water 


vapor, and methane and on monthly averaged climato- 
logical values for other radiatively active trace species. 
These monthly stream function estimates are interpo- 
lated to daily values, which are used to run 2-D back 
trajectories to the 100-hPa surface. Trajectory end- 
points were spaced on a nonuniform vertical grid about 
1 km apart (whose levels are those used in the radia- 
tive calculation; note tick marks on the ordinate) and 2° 
apart in latitude, from 12°S to 12° N. When a trajectory 
crosses 100 hPa, a mixing ratio is assigned to it using 
the saturation mixing ratio given by the minimum zon- 
ally varying temperature at the crossing latitude within 
a week of the crossing date, using 100-hPa temperatures 
from the UKMO. The minimum 100-hPa temperature 
on the latitude circle was chosen in anticipation of the 
likelihood that the coldest tropopause regions, wherever 
they might occur, would be the entry points for water 
vapor into the stratosphere. The mixing ratio results 
thus obtained for 1992 and 1993 were then averaged to 
reduce interannual variablity, particularly that associ- 
ated with the QBO. 

The smooth picture of the time-varying tape signal 
calculated this way (Plate Id) bears a strong resem- 
blance to the satellite observations (Plates la to lc). 
As was shown by Mote et al. [1995], MLS minima and 
maxima at 46 and 22 hPa occur within a month or two 
of the minima and maxima diagnosed by the TEM 2-D 
trajectory calculation. The TEM calculation indicates 
transit times from the tropopause to 46, 22, and 10 hPa 
of 8.5, 13.5, and 17 months, respectively, plus or minus 
about 2 months. 

For a clearer picture of the vertical structure im- 
plied by Plate 1, we show in Figure 1 profiles of q 
from HALOE and from the NOAA Lyman a hygrome- 
ter aboard the NASA ER-2 research aircraft in north- 
ern hemisphere autumn 1994. Data from both upward 
and downward paths of the ER-2 are show T n, the scat- 
ter indicating the spatial variability at a pressure level. 
The HALOE profile represents an average of 7 tropi- 
cal profiles in early November 1994. There are several 
differences between the aircraft and the satellite mea- 
surements that are consistent with the data character- 
istics. The HALOE instrument has approximately a 2 
km resolution (the data points, at 0.2 km intervals, rep- 
resent layers 2 km thick), while the NOAA hygrometer 
has a much finer vertical resolution. In addition, these 
altitudes are near the lower limit for HALOE! q , and 
some questionable data may be included in the aver- 
age. HALOE water vapor tends to be generally lower 
than Lyman a water vapor, a bias that is not fully un- 
derstood. 

Both HALOE and the Lyman a hygrometer show 
a minimum at about 100 hPa, a maximum at about 
70 hPa in air that presumably left the tropopause a 
few months earlier, and a decrease above toward a sec- 
ondary minimum in air that presumably left the tropo- 
pause the previous northern winter. Corresponding 
HALOE and ER-2 profiles in other seasons (not shown) 
have a vertical structure similarly consistent with the 
tape recorder hypothesis; for example, in July, q values 
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Figure 2. Comparison at (a) 22 hPa and (b) 46 hPa daily tropical (I0°S to 10° N) q measurements 
for version 3 MLS (solid curve), for new nonlinear MLS (+ ), for HALOE (A), and for CLAES 
(dashed curve) during the period October 1991 to May 1995. Tick marks are at the beginnings 
of the months indicated. 


at the tropopause are closer to 4 ppmv, the profile mini- 
mum occurs at about 65 hPa, and a maximum occurs at 
about 35 hPa. Further corroboration comes from bal- 
loon profiles (J. Ovarlez, personal communication, 1995; 
Ovarlez et ah, submitted manuscript, 1995) in January- 
February 1994, when the profile maximum occurred at 
about 55 hPa. Weinstock et al. [1995] compared aircraft 
profiles taken in different months by using a temporally 
constant diabatic vertical velocity to “rewind the tape” 
to the previous January 1, as we would say it in terms of 
the tape-recorder analogy. The results from both these 
studies corroborate the tape-recorder hypothesis. 

In Figure 2 we compare daily mean tropical mea- 
surements of q by MLS, HALOE, and CLAES, at two 
lower stratospheric levels for nearly 4 years after the 
launch of UARS. We show both the standard version 3 
for MLS, the solid curves, and the results of the new, 
nonlinear retrieval, the plus symbols. For the period 
when all three instruments were operating, their an- 
nual variations are in moderately good agreement at 
22 hPa (Figure 2a). HALOE q continues the pattern, 
with the lowest values in about April and the highest 
in about November, though with large interannual am- 
plitude fluctuations. CLAES and MLS exhibit periodic 


irregularities associated with the yaw maneuver that 
UARS must perform. There are systematic differences 
among the instruments, however. Version 3 MLS water 
vapor is uniformly about 0.5 ppmv higher than HALOE 
and CLAES; this is a known bias of MLS in the middle 
and upper stratosphere, and work is under way to cor- 
rect it. The nonlinear retrieval produces values that are 
0. 1-0.2 ppmv lower than those from version 3, but still 
higher than those from HALOE and CLAES. HALOE 
q values are initially much lower than CLAES q but 
converge with time toward CLAES q. 

At 46 hPa (Figure 2b), the instruments agree less 
well, especially in 1992. CLAES and HALOE may be 
significantly affected by the presence of volcanic aerosol, 
which decreases with time; CLAES data before July 
1992 are omitted, and HALOE data before April 1992. 
The annual cycle in HALOE q is somewhat clearer than 
at 22 hPa, with a minimum in about October or Novem- 
ber and a maximum in about April or May. MLS and 
HALOE show similar annual variations (though both 
MLS retrievals are again higher), but MLS version 3 and 
the nonlinear retrieval have somewhat different phases 
of the annual cycle. We believe that this phase shift 
may be partly due to the influence of the annual cy- 
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Figure 3. Seasonal variation of residual vertical ve- 
locity at 78 hPa from October 1991 to October 1994, 
from the TEM calculation described in section 3. Con- 
tour interval is 0.05 x 10" 3 m s" 1 , negative contours are 
dashed, and values greater than 0.3 x 10“ 3 in s -1 are 
shaded. 

cle in the climatology used in the standard version 3, 
which at 46 hPa is in phase with the 100-hPa annual 
cycle. Since the nonlinear retrieval not only seems to 
agree better with IIALOE, but also relies less on clima- 
tology at this level than did the version 3 retrieval (in 
particular using a constant climatological profile rather 
than imposing an annual cycle), we are inclined to place 
more confidence in the results of the nonlinear retrieval. 

3.1. Variations in Ascent Rate 

We noted above that Plate lb implied a larger mean 
ascent rate during northern winter (about 0.4 x 10“ 3 m 
s~ 1 ) than during the rest of the year (about 0.2 x 10“ 3 m 
s _l ). This is also clear in the original TEM 2-D trajec- 
tory calculation (Plate Id). For a more detailed view, 
in Figure 3 we show the latitude-time variation of the 
vertical velocity at 78 hPa, just above the tropopause, 
as calculated by Rosenlof[ 1995]. The annual maximum 
(about 0.5 x 10~ 3 m s“ l ) occurs south of the equator 
in December or January, and large ascent rates (greater 
than 0.3 x 10“ 3 m s” 1 ) also occur north of 10° N between 
March and July. There is also interannual variation of 
the latitude of the southern maximum, which in 1991- 
1992 and 1993-1994 lies further off the equator than in 
1992-1993. 

Carr et ai [1995] found (but did not show) an in- 
triguing correlation between MLS q at 22 hPa and an 
index of QBO-induced anomalies in vertical velocity. 
The index was defined as the zonal wind at 10 hPa 
minus the zonal wind at 50 hPa. Successive minima 
in q at both 22 and 46 hPa appear to be at least 14 
months apart in the MLS record (see Figure 2), raising 
questions about the relative importance of the QBO 
and the tropopause annual cycle in determining lower 
stratospheric q values. In order to clarify the role of the 
QBO we used a two-dimensional model, which wc now 
describe. 


The “THIN AIR” model is described by Kinnersley 
[1995, 1996] and is an improved version of that de- 
scribed by Kinnersley and Harwood [1993]. Its improve- 
ments lie chiefly in the addition of three mutually inter- 
acting planetary-scale (zonal wavenumber T-3) Ilossby 
waves, whose breaking is parameterized and which are 
used to calculate consistently the Ilossby wave drag on 
the zonal wind (important for the wave-driven pump- 
ing of the global circulation) and to parameterize the 
Ilossby- wave- induced flux of trace gases. The wave- 
breaking parameterization is similar but not identical 
to that of Garcia [1991]. The model produces realistic 
zonal winds, temperatures, semiannual oscillation, an- 
nual cycles, and Ross by wave amplitudes [Kinnersley, 
1996]. The parameterized QBO, with a period of al- 
most exactly 2 years, is driven by the wave drag due to 
a thermally and frictionaliy damped Kelvin wave with 
a phase velocity of 20 m s” 1 , and an “anti- Kelvin wave” 
with a phase speed of -20 m s -1 and the same disper- 
sion relation and wave-induced momentum flux as the 
Kelvin wave [Kinnersley and Pawson, 1996]. The model 
is formulated with potential temperature as the verti- 
cal coordinate, except in the troposphere, where hybrid 
coordinates are used to avoid intersecting the ground. 
Levels are spaced at about 3.5 km from the ground to 
100 km. The latitudinal resolution is about 9.5°. 

The model uses an upwind advection scheme, which 
has numerical diffusion that is unsatisfactory for han- 
dling strong gradients; we found that when we intro- 
duced an annually varying boundary condition for q at 
the tropical tropopause, the annual variations traveled 
upward too quickly in comparison with the calculations 
of Mote et ai [1995], even though the vertical veloci- 
ties were similar. We therefore represented the annual 
variations of q by trajectories, which were calculated us- 
ing modeled vertical velocities and ignoring horizontal 
motion except at the model tropopause. Parcels were 
assumed to leave the tropopause from the latitude box 
that had the largest upward velocity, thus accounting 
for the seasonal movement of upwelling (9.5°S in Jan- 
uary to 9.5° N in July; compare with Figure 3), but after 
leaving this box they were assumed to remain over the 
equator and so were advected with the equatorial ver- 
tical velocity. Each parcel was assigned an annually 
varying initial value of q with a minimum of 2 ppmv 
in February and a maximum of 4 ppmv in August, and 
the values of q were relaxed to 3 ppmv with an e-folding 
time of 200 days to represent the attenuation of the an- 
nual cycle with height. It should be noted that the at- 
tenuation time in the real tropics appears to be longer; 
see section 3.3. 

The only source of interannual variability in the THIN 
AIR model is the QBO; and in a 7-year run we found 
alternate years to be nearly identical. Figure 4a shows 4 
years of results from the THIN AIR trajectory calcula- 
tion. The vertical derivative of zonal mean zonal wind is 
also shown (dark contours) and will be discussed below. 
Note that the ascent rate of the water vapor contours 
varies from year to year depending on the sign of the 
wind shear, an effect that is most pronounced between 
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Figure 4. (a) THIN AIR 2-D model equatorial parcel water vapor mixing ratio (light shading, 
moist; dark shading, dry; contour interval, 0.2 ppmv). Also plotted is the vertical shear of the 
zonal mean zonal wind (dashed contours, negative shear; (solid contours, zero or positive shear; 
contour interval, 5 m s -1 km -1 ), (b) Time series of parcel water vapor at 22.5 km (about 40 
hPa) from the THIN AIR model. Intervals between successive minima and maxima (labelled by 
triangles and squares, respectively) are indicated. 


18 and 22 km. Figure 4b shows a slice through Figure 
4a, a time series of mixing ratio of parcels passing the 
22.5 km (about 40 hPa) level. The times between suc- 
cessive minima or successive maxima are labelled; the 
minima are from 10.2 to 13.6 months apart depending 
on the phase of the model QBO. 

Easterly or negative vertical shear of zonal wind 
(dashed contours) implies lower equatorial tempera- 
tures and enhanced ascent, while westerly or positive 
vertical shear implies higher temperatures and reduced 
ascent. These effects are clear in Figure 4a: regions 
of positive shear show slower ascent (flatter contours) 
than regions of negative shear. 

In the real atmosphere from November 1991 to about, 
May 1992, negative shear of zonal wind prevailed be- 
tween 20 and 100 hPa [Siombank and O'Neill , 1994], 
implying earlier arrival of minima and maxima in q at a 
given pressure level. As the QBO easterlies descended, 
giving way to westerlies, the Q BO-induced vertical mo- 
tion caused reduced ascent rates below 22 hPa, implying 
later arrival of minima and maxima in q at 22 hPa. A 1- 
month-early arrival of the 1991 100-hPa minimum, due 
at 22 hPa in March 1992, and a 1-month-late arrival of 
the 1992 100-hPa minimum, due at 22 hPa in March 
1993, explain the 14-month interval between minima 
in the left half of Figure 2a as indicated by MLS and 
CLAES. The model QBO (Figure 4) produces a very 


similar interval (13.6 months) and shows that shorter 
intervals should alternate with longer intervals to pro- 
duce a long-term mean interval of about 12 months. 
The timing is harder to deduce from the longer HALOE 
record because of its poorer time resolution, but the 
1993 100-hPa minimum seems to arrive at 22 hPa in 
about March 1994. The shapes of HALOE vertical pro- 
files (not shown) differ from year to year in a manner 
that is consistent with stretching or compressing of the 
tape signal by the QBO. Note that because the TEM 
calculation was based on an average of 1992 and 1993 
data, there is no interannual variability in Plate Id as 
there is in Figure 4a. 

3.2. Variations in Tropopause Temperature 

To check the tape recorder hypothesis quantitatively, 
we would have to find a way to relate the observed tape 
signal to the mixing ratio of air entering the strato- 
sphere. In formulating our TEM 2-D trajectory calcu- 
lation, we assumed that the effective entry mixing ratio 
could be found by tracing tropical air back to the 
latitude at which it crossed the 100-hPa surface and 
finding the saturation mixing ratio q s not at the mean, 
but at the lowest UKMO 100-hPa temperature on the 
corresponding latitude circle. We now examine several 
facets of this calculation, and the detailed assumptions 
involved, in an attempt to explain the amplitude dis- 
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crepancy mentioned in section 1: 

( 1 ) the selection of an appropriate area of the tropics 
for the air to enter the stratosphere, 

(2) the difference between 100-hPa and profile min- 
imum <? s , 

(3) the difference between gridded and actual tem- 
peratures, 

(4) relevant aspects of the UKMO data assimilation 
and interpolation to UARS levels, 

(5) large-scale 3-D motions not accounted for by the 
TEM calculation, 

(6) small-scale, unresolved 3-D motions near the tropo- 
pause, including the possibility of “hyperventilation.” 
These points are discussed in order in the following sub- 
sections. 

3.2.1. Selection of an appropriate area of the 
tropics for the air to enter. It is generally accepted 
that tropical area- averaged tropopause temperatures 
are too high to explain observed q values [Newell and 
Gould’Stewari , 1981; Holton , 1984]. We must there- 
fore choose our g s values from only some portion of the 
tropics where temperatures are below the mean. The 
first and simplest approach would be to use the lowest 
100-hPa temperature anywhere in the tropics. But one 
might expect this to underestimate q e severely, and in 
any case it seems unrealistic in that it implies an as- 
sumption that all air enters at one grid point. A sec- 
ond approach would be to assume that air enters the 
stratosphere only at times and locations whose satura- 
tion mixing ratio q s falls below some upper limit; this is 
the "‘stratospheric fountain” hypothesis of Newell and 
Gould-Stewart [1981]. A third approach is the one we 
chose: to rely on the TEM calculation to select (and 
apply appropriate weights to) the latitudes where air in 
the tropics first crossed the 100-hPa surface, and to use 
the minimum q s on each latitude circle. This approach 
too has limitations, and might also be expected to un- 
derestimate q E \ hut, we find that the situation is not so 
simple. 

To see what result might be expected from the third 
approach, we consider regional averages of q s . Aiitcks 
and Robinson [1983], using radiosonde data from 1979, 
presented statistics of tropopause and 100-hPa q s in dif- 
ferent periods and regions in the tropics. They showed 
that very low tropopause values of q s , with extremes 
below 1.6 ppmv, were most common during northern 
winter in the region 10°S to 10°N, 90°E to 150°W (ap- 
proximately the “fountain” region of Newell and Gould- 
Stewart [1981]). For this whole region and for the period 
January 25 to February 12, 1979, we find, from Atticks 
and Robinson’s tables, a regional average tropopause 
q s of 3.5 ppmv. The TEM 2-D trajectory calculation 
yields a minimum q E of 2.9 ppmv, lying between the re- 
gional average 3.5 ppmv and the extreme values below 
1.6 ppmv and remarkably close, moreover, to the min- 
imum q values observed by HALOE at 46 hPa, around 
3.0 ppmv (Figure 2b). But in northern summer the 
situation is completely different. For the period July 
5 22, 1979, the average q s for the same region is 6.1 
ppmv. The TEM calculation yields a maximum q E of 


7.5 ppmv. Both values are well above the maximum 
observed at 46 hPa, around 4.2 ppmv (Figure 2b). It 
appears that the TEM calculation is grossly overesti- 
mating, not underestimating, q e in northern summer. 

In addition, Atticks and Robinson’s tables allow us to 
calculate the regional average tropopause q s for July 5- 
22 in the Asian monsoon region (10°-25°N, 70-1 10°E). 
This regional average is found to be 4.7 ppmv, closer 
to the 46 hPa maximum of about 4.2 ppmv. Now, in 
the TEM calculation the 2-D back trajectories from the 
middle stratosphere between 10°S and 10° N generally 
cross the tropopause in a narrower latitude band, usu- 
ally south of 12° N. In the TEM zonal mean picture, 
taken literally, air from the colder tropopause over the 
Asian monsoon region passes upward, northward, then 
downward again, such that it never reaches the mid- 
dle stratosphere (see the plot of residual mean stream 
function in July presented by Roseniof [1995]) and thus 
does not influence the TEM calculation of q E . If this 
were truly what happens to air parcels (a point to which 
we return in sections 3.2.5 and 3.2.6), then the region 
of the monsoon would have no influence on middle- 
stratospheric trace constituents, including water vapor, 
between 12°S and 12°N. 

However, even if the TEM trajectories were replaced 
by more realistic, 3-D trajectories, there would still be 
another problem. Suppose for instance that all the air 
at 46 hPa that passed through the tropical tropopause 
in northern summer came from the Asian monsoon re- 
gion. The 100-hPa UKMO temperatures in that re- 
gion would be inconsistent with the notion that the air 
was dried to the saturation mixing ratio q s there. The 
100-hPa temperatures for which q s equals the lowest 
(3.0 ppmv) and highest (4.2 ppmv) reliable values of q 
observed at 46 hPa by UARS instruments (Figure 2b) 
are 189. 6K and 191. 6K, respectively. These values are 
rather rare in UKMO 100-hPa temperature data. This 
can be seen from Figure 5. Figure 5a shows time se- 
ries of the average and minimum UKMO 100-hPa tem- 
peratures between 30°S and 30°N (note the horizontal 
line showing 189.6 K). Figure 5b shows the fraction of 
the area of the 30° S to 30° N region with temperatures 
less than 191.6 K. Assuming that the driest air at 46 
hPa entered the stratosphere in northern winter, only 
a few grid points could explain a value of 3.0 ppmv, let 
alone lower values. (In fact both HALOE and the MLS 
nonlinear retrieval indicate much lower values than this 
at levels between 100 and 46 hPa, though the signifi- 
cance is unclear because the data there have larger error 
bars.) Also, although some grid points (generally fewer 
than 5% of them) are cold enough in northern winter 
to explain a value of 4.2 ppmv, none are cold enough in 
northern summer, when the minimum 100-hPa q s ex- 
ceeds 7 ppmv. It is clear from Figure 5, then, that any 
air entering the stratosphere in northern summer and 
dried to saturation (at the UKMO 100-hPa tempera- 
tures) would be much moister than observed. 

Is it possible that air simply does not enter the strato- 
sphere when q s is too high, that is, for several months 
in northern summer? On the contrary, the C0 2 data of 
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Figure 5. Statistics from 100-hPa UKMO assimilated 
temperature data for October 1991 to April 1994, for 
the area between 30° S and 30° N. (a) Area average tem- 
perature (top curve) and minimum temperature (bot- 
tom curve); horizontal line indicates 189.6 K. (b) Frac- 
tion of the area covered by temperatures less than 191.6 
K. The values 189.6 K and 191.6 K are the temperatures 
at which q s at 100 hPa equals the lowest and highest 
reliable q values in Figure 2, taken respectively as 3.0 
ppmv and 4.2 ppmv. 

Boering et al. [1995] suggest continuous entry of air to 
the stratosphere, a view consistent with the satellite ob- 
servations of q (Plate 1). Furthermore, the TEM calcu- 
lation has nonzero vertical velocity at 100 hPa through- 
out the year, as also expected for robust dynamical rea- 
sons [e.g., Holion et al., 1995, Haynes et al, 1996]. It 
therefore seems more sensible to allow the TEM cal- 
culation, as far as it is able, to select the appropriate 
region, rather than to use only 100-hPa temperatures 
that fall below a certain limit, which in any case would 
not make sense in northern summer because then the 
limit is not attained at all. Our choice of the minimum 
q s on each latitude circle, though apparently somewhat 
arbitrary, did succeed remarkably well in compensating 
the errors associated with points 2, 3 and 4 in north- 
ern winter though not in northern summer. One factor 
to which the northern summer discrepancy points very 
strongly is the difference between 100-hPa q s and profile 
minimum q s . This is considered in the next subsection. 

First, however, we should ask how important inter- 
annual variability might be. Is it reasonable to com- 
pare area-averaged q s values for 1979 with data from 
different years, in particular the qv values diagnosed 
from 1992-1993 100-hPa temperatures? Frederick and 
Douglass [1983] showed annual mean 100-hPa temper- 
atures for 1972-1979 for several tropical stations. For 
most such stations, especially those with lowest mean 
temperatures, interannual variations were of the order 
of 1 K or less. They also stated [Frederick and Dou- 
glass, 1983, p. 1399] that “most of the variance [about 


the long-term mean] is in the annual cycle.” Our own 
inspection of 100-hPa analyses from the European Cen- 
tre for Medium- Range Weather Forecasts (ECMWF), 
for 1986-1989 suggests that some interannual variation 
of 100-hPa temperatures in the western Pacific results 
merely from the movement of the temperature contours, 
which has little effect on the area-averaged <? s . Figure 
5 also suggests that interannual variations are smaller 
than annual and intra-annual variations. 

3.2.2. Difference between 100-hPa q s and pro- 
file minimum g s . Frederick and Douglass [1983] have 
noted that if the tropopause pressure p t is less than 100 
hPa, then replacing tropopause temperature T t by 100- 
hPa temperature 7\oo would overestimate profile mini- 
mum q S) because the higher temperature gives higher q s , 
whereas replacing p t by 100 hPa would underestimate 
r/ s , because q s oc p -1 where p is air pressure. It hap- 
pens that the two effects approximately cancel when 
pt < 100 hPa. Conversely, the two effects reinforce 
when the p t > 100 hPa, as is common in northern sum- 
mer, giving an overestimate of minimum q s . Retd and 
Gage [1981] showed that, over a large area of the trop- 
ics, tropopause heights are indeed lowest in northern 
summer and highest in northern winter. The question, 
then, is: could the seasonal migration of the tropopause 
across the 100-hPa surface explain the amplitude dis- 
crepancy in the tape signal? 

We look first at some relevant observational infor- 
mation for the equatorial western Pacific, yielding time 
series of the two types of q s estimate just discussed, say 
qf sloo and q st , q s 100 being the estimate from Tioo and 
q st being the estimate from Tt and pt- Time series of 
monthly mean 7ioo> T t , and p t derived from radiosonde 
profiles were kindly provided by G. Reid for the island 
of Koror in the equatorial western Pacific (7°N, 134°E), 
for the period 1965-1981. Figure 6a compares the re- 
sulting <7sioo an< l ?st> an d Figure 6b shows the difference 
(?sioo “ 7st)* From November to May the tropopause is 
between 95 and 98 hPa and the mean thermal structure 
is such that q s \oo w 7s t ■ This illustrates the approx- 
imate cancellation already mentioned. From June to 
September, however, the difference between tfsioo auc ^ 
q st is much larger, close to 2 ppmv. Atticks and Robin- 
son [1983] display an area mean vertical profile of tem- 
perature for the 150° E to 150° W sector, just east of 
Koror, and there too qr sl00 & q st in northern winter but 
not in northern summer. 

It is possible of course that the use of time series 
of monthly means may obscure some relevant informa- 
tion. Also, as Frederick and Douglass noted, the p- 
dependence of q s implies that there is an additional dif- 
ference, which is usually small, between the q st defined 
above and the true profile minimum value of q s . 

Figure 6 goes some way (2 ppmv) toward explain- 
ing the amplitude discrepancy in the TEM 2-D tra- 
jectory calculation, but are such differences between 
q s \oo an d Qst found throughout the tropics? In north- 
ern winter, the lowest T t values usually occur near 
Koror, but in northern summer they occur over the 
Asian monsoon region, and a number of studies show 
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Figure 6. Seasonal variation of the error from using 
100-hPa temperatures to estimate minimum q s , illus- 
trated by monthly mean radiosonde data from Koror 
(7°N, 134°E) for 1965-1981. (a) Solid curve shows the 
estimate g sloo minimum q s , computed from 100-hPa 
temperature and pressure; dotted curve shows the es- 
timate g 8t from tropopause temperature and pressure, 
(b) Curve shows the difference (q s 100 — q 8t ) on an ex_ 
panded scale. Data by courtesy of G. Reid. 


that p t < 100 hPa there. Atticks and Robinson [1983] 
binned profile-minimum temperature T mm and the cor- 
responding pressure p m j n for hundreds of radiosonde as- 
cents. For the region 10°S to 10°N, 150°E to 150° W 
in northern winter, the bin with the most data points 
had p m in = 92.5 hPa, and for the region 10°-25°N, 
70°-110°E in northern summer, the bin with the most 
points had p mm = 97.5 hPa. Also, Krishna Murthy 
et al. [1986] showed that over Delhi (but not south- 
ern India), the monthly mean p t was less than 100 hPa 
in summer. Finally, an analysis of ECMWF data (E. 
Highwood, personal communication, 1995), which used 
lapse rates to find the tropopause between vertical grid 
points, likewise showed p t less than 100 hPa over much 
of northern India and the Bay of Bengal in northern 
summer. These studies imply that, in the Asian mon- 
soon region during northern summer, g SH)0 c °uld be a 
reasonable proxy for q st if gridded Tioo values are ac- 
curate enough. The 2 ppmv difference during northern 
summer (Figure 6b) probably does not apply to the 
region of lowest T* (the Asian monsoon) then. There 
is, however, an unexplained discrepancy between the q s 
values (> 7 ppmv) diagnosed from UKMO minimum 
100 hPa temperatures (Figure 5a) in northern summer 
and area averaged q s values (4.7 ppmv) derived from 
the tables of Atticks and Robinson [1983]. 

3.2.3. Difference between gridded and actual 
temperatures. 

The caveat about gridded Tioo being accurate enough 
is crucial. For instance, even over the Asian monsoon 


region, the UKMO gridded Tioo values are far too high, 
as mentioned earlier, to account for observed q values 
at 46 hPa: the gridded values gave q s exceeding 7 ppmv 
as compared with the 46 hPa maximum of about 4.2 
ppmv. One question is whether horizontal grid reso- 
lution limitations are important, especially as there is 
ample reason to believe that cumulonimbus convection 
plays an important role, locally, in the process of trans- 
ferring mass to the stratosphere [e.g., Kelly et al , 1993; 
Kritz et al, 1993]. Processes at and near the tops of 
penetrating convective clouds are probably an essen- 
tial component especially as regards estimating ?e [e.g., 
Danielsen, 1982]; but the extremely low temperatures 
involved are both too transient and too fine in scale to 
be well represented even in radiosonde data, let alone 
in gridded data after assimilation. 

Despite these limitations we may try to see what 
radiosonde data can tell us about the frequency and 
seasonal distribution of very low temperatures and q s 
values, and try to relate these to our grid-based esti- 
mates of <?e* The study by Atticks and Robinson [1983] 
is again helpful. They plotted T mi n against p mm for 
each of hundreds of radiosonde profiles. For the region 
10°S to 10°N, 150°E to 150°W over a 20-day period 
in northern winter, dozens of profiles had T min values 
less than 186 K (the lowest UKMO Tioo value in Fig- 
ure 5), and at p m in values ranging from 80 to 110 hPa, 
implying g Smin approaching 1 ppmv. They noted that 
none of these radiosonde profiles coincided with active 
convection, so it is possible that convective events would 
briefly produce even lower q Smm values. It seems plausi- 
ble therefore that the synoptic scale of the gridded data 
is insufficient to resolve the smaller scales, where very 
low temperatures associated with convection are proba- 
bly of key importance to qv- The implication is that, for 
the purpose of estimating q&, a warm bias is inherent 
in the gridding process. In our TEM 2-D calculation 
we sought to compensate, in effect, for the warm bias 
by choosing the lowest available temperature on each 
latitude circle. 

3.2.4. Relevant aspects of the UKMO data as- 
similation and interpolation to UARS levels. As 

the preceding sections suggest, the vertical and hori- 
zontal resolution of assimilated data may sharply limit 
our ability to estimate </e- Bor the TEM 2-D calcu- 
lation we used UKMO data because they were conve- 
niently available as part of the UARS correlative data 
base and they were on a regular spatial grid, simpli- 
fying the calculation of a weighted-average q s as our 
estimate of <je- In the statistics produced as part of 
the UKMO data assimilation process, the region be- 
tween 30° S and 30° N and between 50 and 100 hPa has 
one of the largest absolute means (-~0.5K) of the quan- 
tity defined as observed minus analyzed temperature 
(R. Swinbank, personal communication, 1994). This is 
an additional warm bias. 

Another problem is that Tioo values have to be found 
by interpolating between the 89 and 112 hPa model 
levels. Although the interpolation is nonlinear, it is 
likely to overestimate T mm . The combination of an 
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assimilation bias and an interpolation bias, adding to 
the horizontal-resolution bias discussed in section 3.2.3, 
could well be significant: if the minimum temperature 
curve in Figure 5a were shifted down by, say, a degree, 
the fraction of grid points cold enough to explain ob- 
served mixing ratios would increase markedly. 

The bias in UKMO tropical lower stratospheric tem- 
peratures can have an effect on the TEM stream func- 
tion itself. Uncorrected, it would yield vertical velocities 
that are too weak. But because we correct for global 
mass balance, we can reasonably assume that the worst 
such errors have been taken out, though there remain 
other complications [e . g . , Shine , 1989], and it is difficult 
to be sure just what the net effect would be. 

3.2.5* Large-scale 3-D motion not accounted 
for by the TEM 2-D calculation. Our approach 
to calculating by backward trajectories to the tropo- 
pause is limited by the accuracy of the TEM calculation 
in the lower stratosphere. Radiative errors, though esti- 
mated to be small, arise from many sources besides the 
important temperature errors already discussed. Out- 
side the radiation scheme itself there is the simplified 
climatological treatment of clouds, and the errors in 
measurements of ozone and water vapor in the lower 
stratosphere. It seems plausible that such errors could 
contribute to errors in the source latitude of strato- 
spheric air. But a still more serious error arises from 
zonal averaging. Stratospheric air at a given latitude 
could very well trace its origins to a far wider range of 
latitudes than the TEM 2-D calculation would indicate; 
and this could be especially serious in northern summer, 
as already noted. 

One of us (K.H.R.) conducted a preliminary study of 
the impact of such errors by running 3-D backward tra- 
jectories from 55 hPa using observed (UKMO) winds, 
starting in January and July and running the trajecto- 
ries back for 6 months. These results show that air does 
indeed come from a wider range of latitudes than the 
TEM calculation would indicate, and furthermore that 
it crossed the tropopause at a range of times. Many of 
the parcels that started in January crossed the tropo- 
pause in the Asian monsoon region during the previ- 
ous summer, but many also crossed the tropopause in 
September or October. Even more came from the ex- 
tratropical stratosphere. 

The fact that the (mass balance corrected) TEM cal- 
culation nevertheless produces remarkably good phase 
agreement with observed q (Plate 1; Mote et al. [1995]) 
argues strongly for corresponding accuracy of the TEM 
vertical velocities. Indeed, Plate 1 can be looked on as 
a stringent independent check on such accuracy. What 
the TEM 2-D calculation cannot adequately represent 
is all the detail in the highly disturbed layers just above 
the tropopause, hence statistics of pathways from the 
tropopause, hence q^ values. 

3.2.6. Small-scale, unresolved 3-D motion near 
the tropopause, including the possibility of “hy- 
perventilation.” It is possible that gridded data mis- 
represents not only temperatures but also transport. 
For example, suppose that there exists a sub-grid-scale 


mechanism for selecting only the driest air to be taken 
up into the tropical middle stratosphere. One such 
mechanism that is physically possible, indeed arguably 
likely to operate for part of the time at least, is what 
might be called “hyperventilation” of the tropical lower 
stratosphere by the tallest cumulonimbus clouds. Hy- 
perventilation means that the tallest clouds inject sig- 
nificantly more dry air into the tropical lower strato- 
sphere than can be taken up into the middle strato- 
sphere by the large-scale (extratropically controlled) ad- 
vection. In a hyperventilated lower stratosphere, the ex- 
cess mass injected has to be displaced back downward 
or sideways. 

The possibility of hyperventilation arises from the 
possibility that the tail of the statistical distribution of 
tall cumulonimbus injection altitudes is associated with 
a strongly negative contribution to the subgrid-scale up- 
ward eddy flux of water vapor in the lower stratosphere, 
while not being sufficient to carve out a new tropopause 
at the same altitude. The eddy flux could behave in 
this way because the upward motion associated with 
tall cumulonimbus penetrations is likely to be strongly 
correlated with the lowest extreme q values, i.e. with 
strongly negative q' . The water vapor fluxes can thus 
behave differently from the eddy heat fluxes influencing 
tropopause formation. 

A hyperventilated tropical lower stratosphere would 
evidently imply effective temperature and <?e values 
substantially less than those calculated from a larger- 
scale view. The question of whether and when such 
hyperventilation is significant for values in the real 
atmosphere needs, and deserves, further investigation. 
This will have to await detailed studies of cloud statis- 
tics, particularly the total mass flux from the tallest 
clouds into the stratosphere. 

3.3. Mixing 

The observed tape signal contains useful new infor- 
mation, difficult to obtain otherwise, about vertical 
mixing within the tropics and near-horizontal isentropic 
mixing into the tropics. This complements what is al- 
ready known about mixing out from the tropics from, 
for instance, the studies of nuclear bomb test debris 
[Feely and Spar , 1960] and of water vapor [ Masten - 
brook and Oilmans , 1983; McCormick et aL, 1993] and 
from diagnoses of transport using analyzed winds [Chen 
et ai, 1994; Waugh, 1996]. These studies have indi- 
cated that near-horizontal isentropic transport out of 
the tropics is more rapid and effective in the lowest few 
kilometers of the stratosphere than at somewhat higher 
altitudes. The tape signal provides new information 
about this. 

Mixing-in rates need not be closely comparable to 
mixing-out rates. The very existence of the tape signal 
reminds us that the tropical stratosphere is somewhat 
isolated from the surrounding extratropics [e.g., Trepte 
and Hitchman , 1992], to an extent that varies with alti- 
tude, and that this is presumably due to the existence, 
in the subtropics, of the kind of flexible, quasi-elastic 
“eddy transport barrier” that has been more exten- 
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sively studied at the edges of the winter stratospheric 
polar vortices [e.g., Norton, 1994; Polvani et al. , 1995; 
Holton et ai, 1995]. Such eddy transport barriers de- 
pend on monotonic isentropic gradients of potential vor- 
ticity and horizontal shear acting together in a subtle 
way [Juckes and McIntyre , 1987], and there is no rea- 
son why mixing in to the tropical overworld need be 
as strong as mixing out to the extratropical overworld. 
Both tend to be inhibited by the eddy transport barrier, 
but not necessarily to the same extent because eddy de- 
tails tend to be different (e.g., mixing out is associated 
with "‘erosion” of the barrier by midlatitude eddies). 
Therefore we cannot be certain that rates of mixing in 
will be the same as rates of mixing out. This makes esti- 
mates of mixing in, which are hard to obtain, especially 
interesting. 

In this connection it is useful to note first that H = 
2 (CH 4 )+(H 20 ) values in the extratropical overworld 
are nearly uniform at about the mean value of the trop- 
ical tape signal, about 6.5 ppmv. Note also that extra- 
tropical overworld fluctuations in H are probably due 
mostly to poleward erosion of tropical air and hence 
should be within the range of the tape signal. Larger 
departures from the mean value are found in the meso- 
sphere and in the southern polar vortex, but these are 
unlikely to influence air that reaches the tropics [Mote, 
1995]. Mixing in from the surrounding extratropics, 
therefore, is overwhelmingly likely to attenuate, not am- 
plify, the tape signal. 

Thus motivated, we have attempted to estimate the 
timescales for attenuation of the tape signal using HA- 
LOE H and MLS q (which is better sampled, though 
more complicated to interpret because of methane ox- 
idation). Two methods were used; the first was ap- 
plied both to HALOE H and to MLS q , and the second 
was applied only to HALOE H . It turns out that the 
answers are nearly the same in each case. The first 
method is to calculate a peak-to-peak amplitude using 
the maximum and minimum for the entire data record 
at each altitude. The second, which can be applied only 
to HALOE because of its longer record, is to Fourier 
analyze 3 years of data and take the annual compo- 
nent. The amplitude calculated by the second method 
is smaller than that calculated by the first, but the frac- 
tional attenuation with height, i.e., the ratio of ampli- 
tudes at different heights, turns out to be nearly the 
same. 

In view of the altitude dependence of mixing out al- 
ready mentioned, which presumably indicates altitude 
dependence of the strength of the eddy transport bar- 
rier, we might expect the attenuation of the tape signal 
to show altitude dependence as well. In fact, we found 
that the UARS level at 46 hPa roughly delimits two 
different regimes. Between 46 and 10 hPa the peak-to- 
peak amplitude of both HALOE H and MLS q drops 
from 0.9 ppmv to 0.5 ppmv during the 9 months re- 
quired for the journey. Assuming exponential decay, 
this implies an e-folding timescale r of about 15 18 
months. Between 100 and 46 hPa, however, the e- 
folding timescale r is only about 7-9 months. The 


slower degradation of the tape signal between 46 and 
10 hPa indicates less mixing in from the extratropics 
and/or less vertical mixing. 

The e-folding attenuation times just estimated im- 
ply upper limits to vertical mixing rates in the tropical 
lower stratosphere. There is good reason to suppose 
that vertical mixing in the lower stratosphere can be 
represented in the ensemble mean by a true Fickian 
eddy diffusivity [ Dewan , 1981]. There can be an addi- 
tional (also Fickian) diffusivity from the quasi-random 
fluctuations in diabatic heating experienced by an air 
parcel. If L is a typical vertical length scale, then 
Fickian diffusion implies a diffusive timescale Tt > iff ~ 
K z /L\ 

Substitution of the observed values of r into the re- 
lation tdiff ~ K z f L 2 gives upper limits on the total 
K z from all such effects. This is because such substi- 
tution amounts to assuming zero mixing in from the 
extratropics. Using the data from Plates la and lb to 
estimate L, we get vertical wavelengths (= 2 7r x L) of 
about 5-7 km between 100 and 46 hPa, and 8-12 km 
between 46 and 10 hPa. If there were no mixing-in, 
then we would have ~ t ~ 7-9 months between 

100 and 46 hPa, and this gives a K z about 0.03-0.07 
m 2 s _1 and a best estimate of about 0.04 m 2 s -1 . With 
fDiFF ~ t ~ 15-18 months between 46 and 10 hPa, this 
gives K z about 0.03-0.09 m 2 s _1 and a best estimate of 
about 0.06 m 2 s _1 . 

The independent studies of aircraft data by Avallone 
and Prather [1996], Boering et al [1995], and K. Min- 
schwaner et al (The bulk properties of isentropic mixing 
into the tropics in the lower stratosphere, submitted to 
Journal of Geophysical Research , 1995) indicate that 
some extratropical air is mixed into the tropics, with 
timescales similar to those we have just estimated. If 
so, the actual K z values must be even smaller. Even 
the K z values just deduced, which are upper limits on 
the order of magnitude of K z , as just emphasized, seem 
distinctly smaller than values for K z like 0.2 m 2 s _l es- 
timated from radar studies and from model studies like 
that of Dewan [1981] (also World Meteorological Orga- 
nization (WMO) [1986, and references therein]). 

A further reason for supposing that I\ z values are 
still smaller than the above estimates is the altitude de- 
pendence of the tape signal attenuation. The increase 
in attenuation times above the first few kilometers of 
the tropical stratosphere seems to mirror the upward 
decrease in poleward erosion rates found by Chen et ai 
[1994] and Waugh [1996] from the lower to the middle 
subtropical stratosphere. Even though, as was already 
mentioned, the eddy transport barrier need not inhibit 
mixing in to the same degree as mixing out, all the re- 
sults under discussion seem consistent with hypothesiz- 
ing a more effective subtropical eddy transport barrier 
in the middle stratosphere than in the first few kilome- 
ters of the overworld, i.e. in the first few kilometers 
above the 400 K isentrope. Long 3-D trajectory runs 
by one of us (K.H.R.) were used to trace the origin of 
tropical air at various altitudes, and these results, also, 
indicate that there is an altitude range in the tropics 
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(within the overworld, between about 50 and 20 hPa) 
in which mixing in of midlatitude air is relatively slow. 

4. Conclusions 

Air transported upward by the extratropically driven 
global-scale mean circulation retains, for 18 months 
or more, a distinct “memory” of the tropical tropo- 
pause conditions it encountered. As was first shown 
by Mote et al. [1995], the timings of individual min- 
ima and maxima in tropical stratospheric water vapor 
q at 46 and 22 hPa are reasonably well represented by 
a 2-D transformed Eulerian mean (TEM) back trajec- 
tory calculation. This not only corroborates the tape 
recorder hypothesis but also validates the mass balance 
corrected radiative calculation in the lowest few kilome- 
ters of the tropical stratosphere, where heating rates are 
very small. The longest-running picture of the tape sig- 
nal is provided by HALOE =2(CH 4 )+(H 2 0) (Plate 
lb), which also shows that, in the lowest few kilometers 
of the tropical stratosphere, ascent rates are smaller 
in northern summer (about 0.2 x 10“ 3 m s -1 ) than 
in northern winter (about 0.4 x 10“ 3 m s _1 ). These 
numbers provide a new independent check on the or- 
ders of magnitude of large-scale tropical stratospheric 
upwelling velocities previously estimated in other ways 
[c.g., Holton et al, 1995]. 

From Plate lb it is also clear that, most of the time, 
a tropical profile of q would show two minima and two 
maxima between 100 hPa and 10 hPa, like an upward 
extension of Figure 1. The upward advection of a sea- 
sonal minimum is a necessary part of this picture and 
thus seems to be the most likely explanation of the 
“elevated hygropause” [Kley et al, 1982; Kelly et al, 
1993] than zonal advection from a contemporary distant 
source. 

However, it is not clear in full detail how the initial 
marks on the tape are made, i.e., how the real atmo- 
sphere determines the effective “entry mixing ratio” q& 
of air entering the tropical stratosphere. It became clear 
that establishing a quantitative link between qE and as- 
similated temperatures near the tropopause would be 
an extremely formidable problem, beyond the scope of 
this study if not beyond present capabilities altogether. 
One would need accurate knowledge of many 3-D tra- 
jectories and of sub-grid-scale processes near the tropo- 
pause, including cloud-physics processes, as recalled be- 
low. 

There is a large discrepancy between our q& estimates 
from the TEM 2-D back trajectory calculation and the 
observed mixing ratios q for air at higher altitudes such 
as 46 hPa, which probably entered the stratosphere in 
northern summer. We noted several possible reasons 
for this: 

1. The TEM 2-D back trajectory calculation may 
not be able to trace trajectories back to the right loca- 
tions near the tropopause, and almost certainly fails to 
do so in northern summer. Because the TEM calcula- 
tion treats the tropics as zonally symmetric, it cannot 


describe zonally asymmetric circulations like the an- 
ticyclonic outflow from monsoon regions, which seem 
likely to be important to the determination of q e. In 
the tropical middle stratosphere the TEM calculation 
shows, probably unrealistically, no air originating from 
the Asian monsoon region during northern summer. In 
that region at that time, area- averaged tropopause q s 
(from radiosonde data) is about 4.7 ppmv, close to the 
maximum q observed at 46 hPa, about 4.2 ppmv. Com- 
pounding the problem with the trajectories is the prob- 
lem discussed in section 3.2.2 (see also Figure 6), that 
tropopause q s , hence can be grossly overestimated 
when tropopause pressures exceed 100 hPa, as may hap- 
pen at tropopause locations indicated by unrealistic tra- 
jectories. This appears to account for a substantial part 
of the amplitude discrepancy, that is, the discrepancy 
in the 2-D estimates of qE- 

2. The use of UKMO gridded 100-hPa temperatures 
to estimate qE has shortcomings. In the data assimila- 
tion procedure, information about subgrid-scale tem- 
perature variations (which may be of critical impor- 
tance to the question of how air enters the stratosphere) 
can be lost. Furthermore, the interpolation from model 
levels to 100 hPa is almost certain to overestimate the 
model's tropopause temperature. 

3. Air entering the stratosphere in northern summer 
is more likely to be modified by subsequent convection 
events, owing to the lower tropopause and slow mean as- 
cent rate in the tropical lower stratosphere during that 
time of year. Such “hyperventilation” would replace 
moister air by drier air, like shuffling a horizontal deck 
of cards, leaving a drier signal on the tape than would 
be indicated by the synoptic-scale temperatures. 

4. Moist air entering in northern summer may also be 
diluted by drier air from the extratropics, where mixing 
ratios are then lower. There is evidence from aircraft 
data that some air is mixed in from the extratropics 
between about 100 and 50 hPa, where Plates la and lb 
both show the greatest reduction in tape signal ampli- 
tude. 

Establishing the relative importance of these factors 
would require considerable further work. Although a 
better estimate of the profile minimum q s could be 
obtained by examining several hundred radiosonde as- 
cents, as was done by Atticks and Robinson [1983], the 
improved vertical resolution would come at the expense 
of the horizontal grid and dynamical balance provided 
by a data assimilation system. Furthermore, it is not 
even clear how radiosonde temperatures are related to 
the actual water vapor mixing ratio of air entering the 
stratosphere, since that air may initially contain some 
ice crystals that may evaporate (thereby moistening the 
air), grow by turbulent mixing in a radiatively desta- 
bilized cloud [ Danielsen , 1982], or mostly precipitate 
before evaporating [Knollenberg et al, 1993]. 

HALOE H and MLS q data shed new light on the 
role of mixing in modifying the tape signal. In partic- 
ular, our analysis implies that vertical mixing is weak 
and that the subtropical stratospheric “transport barri- 
ers” [e.g., McIntyre, 1990, Norton , 1994, Polvani et al, 
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1995] are effective at inhibiting transport into the trop- 
ics. The H and q data demonstrate such inhibition (sec- 
tion 3.3) for the tropical lower stratosphere, especially 
between about 50 and 20 hPa. An important corol- 
lary is the greater relevance, to the modeling of long- 
lived tracer transport and tracer correlations, of “trop- 
ical pipe models” [Plumb, 1996], as contrasted with 
one-dimensional “global-mixing models” [ Holton , 1986, 
M ah bn an et al., 1986, Plumb and Ko , 1992]. Hence 
the situation is unlike the one-dimensional, globally in- 
tegrated, vertically diffusive situation studied most re- 
cently by Plumb and Ko [1992], which arises in the limit 
of fast meridional mixing. Piurnb [1996] suggested on 
the basis of earlier evidence that the real tropical strato- 
sphere is for chemical purposes closer to the limit of slow 
mixing, and he called it a “pipe”; we might add that it 
is rather like a “leaky pipe.” Our estimates of an up- 
per bound for K z and a lower bound for the mixing-in 
timescale (about 0.04 m 2 s^and 7-9 months, respec- 
tively, in the altitude range 100-46 hPa, and about 0.06 
m 2 s“ 1 and 15—18 months in the altitude range 46-10 
hPa) imply that the actual K z is much smaller and/or 
that the mixing-in timescale is much longer. 

As is demonstrated by the THIN AIR model, the 
QBO can change the speed of the tape slightly, re- 
ducing or increasing parcel transit times to the mid- 
dle stratosphere by about a month. Given the zonal 
mean winds during the UARS period, this reduction 
or increase partly explains the small discrepancies be- 
tween the TEM calculation and the timing of minima 
and maxima in the UARS record, particularly MLS. 

Taken together and accounting for the apparent bias 
in MLS measurements, the UARS instruments imply 
an annually varying q% whose annual average is about 
3.5 ppmv, assuming negligible methane oxidation be- 
tween the tropopause and 46 hPa. A value of 3.5 ppmv 
is in agreement with an estimate of 3.25 ppmv using 
LIMS and SAMS data [Hansen and Robinson , 1989] and 
with a careful estimate (2.9 -3.4 ppmv) using HALOE 
data (E. Remsberg et al., Estimates of the water vapor 
budget of the stratosphere from UARS HALOE data, 
submitted to Journal of Geophysical Research , 1995) 
but somewhat lower than recent estimates using air- 
craft data: Dessler et al. [1994] derived a value of 4.2 
ppmv, and Weinstock et ai [1995] derived a value of 
4.4 ppmv. Remsberg et al., however, note that because 
of its history, air at aircraft flight level may not yield 
an accurate estimate of annual mean </e- 
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